Quantum Anomalous Hall Effect in Single-layer and Bilayer Graphene 
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The quantum anomalous Hall effect can occur in single and few layer graphene systems that have 
both exchange fields and spin-orbit coupling. In this paper, we present a study of the quantum 
anomalous Hall effect in single-layer and gated bilayer graphene systems with Rashba spin-orbit 
coupling. We compute Berry curvatures at each valley point and find that for single-layer graphene 
the Hall conductivity is quantized at a xy = 2e 2 //i, with each valley contributing a unit conductance 
and a corresponding chiral edge state. In bilayer graphene, we find that the quantized anomalous 
Hall conductivity is twice that of the single-layer case when the gate voltage U is smaller than 
the exchange field M, and zero otherwise. Although the Chern number vanishes when U > M, 
the system still exhibits a quantized valley Hall effect, with the edge states in opposite valleys 
propagating in opposite directions. The possibility of tuning between different topological states 
with an external gate voltage suggests possible graphene-based spintronics applications. 

PACS numbers: 72.80.Vp,73.43.-f,73.63.-b,72.20.-i 



I. INTRODUCTION 

Graphene is a two-dimensional material with a single- 
layer honeycomb lattice of carbon atoms. Its isolation in 
the past decade has generated substantial theoretical and 
experimental research activity pQ . Experimental fabrica- 
tion methods continue to progress, motivating a recent 
focus on the possibility of utilizing graphene as a mate- 
rial for nanoelectronics. At the same time, spintronics 
has also progressed in recent years. The spin degrees of 
freedom can be manipulated to encode information, al- 
lowing fast information processing and immense storage 
capacity [2 . This article is motivated by recent research 
progress which has advanced the prospects for spintronics 
in graphene. 

A key element of spintronics is the presence of spin- 
orbit coupling which allows the spin degrees of freedom 
to be controlled by electrical means. It was pointed out 
some time ago that the quantum spin Hall effect [3] can 
occur in a single-layer graphene sheets because of its in- 
trinsic spin-orbit (SO) coupling. Instrinsic SO coupling 
induces momentum-space Berry curvatures (which act 
like momentum-space magnetic fields) that have opposite 
sign for opposite spin. However, the intrinsic coupling 
strength was later found to be weak (~ 10 -2 — 10 -3 meV 
[UJ6]) enough to make applications of the effect appear 
impractical. Luckily another type of spin-orbit interac- 
tion known as Rashba SO coupling [7 appears when in- 
version symmetry in the graphene plane is broken. This 
SO coupling mixes different spins, so the spin-component 
perpendicular to the graphene plane is no longer con- 
served. When acting alone, Rashba coupling causes the 
resulting spin eigenstates to be chiral. One appealing fea- 
ture of Rashba SO coupling is its tunability by an applied 
gate field Eq. Unfortunately the field-effect Rashba cou- 
pling strength is also weak ( ~ 10 — 100 jaeV per U/nm 



[4j [6]) at practical field strengths. Recent experiments 
[8j |9] and ab initio calculation [10] have however sug- 
gested that surface deposition of impurity adatoms can 
dramatically enhance the Rashba SO coupling strength 
in graphene to ~ 1 — 10 meV, raising the hope that spin 
transport effects induced by Rashba SO coupling might 
be realizable as experiment progresses. 

The quantum anomalous Hall effect (QAHE) is char- 
acterized by a quantized charge Hall conductance in an 
insulating state. Unlike the quantum Hall effect, which 
arises from Landau level quantization in a strong mag- 
netic field, QAHE is induced by internal magnetization 
and SO coupling. Although there have been a number 
of theoretical studies of this unusual effect [TTHT6] , no 
experimental observation has been reported so far. In a 
recent article [10 , we predicted on the basis of tight- 
binding lattice models and ab initio calculations that 
QAHE can occur in single-layer graphene in the presence 
of both an exchange field and Rashba SO coupling. In 
this paper, we complement our previous numerical work 
with a continuum model study which yields more analyt- 
ical progress and provides clearer insight into our find- 
ings. We also present a more detailed and systematic in- 
vestigation of the topological phases of both single-layer 
and gated bilayer graphene with strong Rashba SO in- 
teractions. In single-layer graphene the Hamiltonian is 
analytically diagonalizable and we obtain an analytic ex- 
pression for the Berry curvature and use it to evaluate 
the Chern number. For bilayer graphene, the possibility 
of a gate field applied across the bilayer introduces a tun- 
able parameter which we show can induce a topological 
phase transition. We find that when the bilayer poten- 
tial difference U is smaller than the exchange field M . the 
system behaves as a quantum anomalous Hall insulator, 
whereas for U > M, the system behaves as a quantum 
valley Hall insulator with zero Chern number. For each 
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case, we also study the edge state properties of the cor- 
responding finite system using a numerical tight-binding 
approach. 

The paper is organized as follows. We first discuss the 
bulk topological properties for the case of single-layer 
graphene in section II. In section III, we turn our at- 
tention to the case of gated bilayer graphene. We then 
discuss the edge state properties of both the single-layer 
and the gated bilayer cases in section IV. Our conclusion 
are presented in section V. An appendix follows that de- 
velops an envelope function analysis of the edge states in 
the single-layer system. 



II. SINGLE-LAYER GRAPHENE 
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The Brillouin zone of graphene is hexagonal with two 
inequivalent K and K' points located at the zone corners. 
The band structure has a linear band crossings at both K 
and K' . At wave vectors near either of these valley points, 
the envelope-function wavefunctions satisfy a massless 
Dirac equation. We represent the graphene envelope- 
function Hamiltonian in the basis {A "f , B "f , A ^, B 1} for 
both valleys K and K'. 

Rashba SO coupling in graphene was first discussed by 
Kane and Mele [3] and subsequently by Rashba [17 , and 
also in a number of recent papers [T8H20] . The Hamil- 
tonian for valley K including Rashba SO coupling and 
exchange field is 



H = vcr ■ kl s + a (a x s) z 



ML 
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where a and s are Pauli matrices that correspond respec- 
tively to the pseudospin (i.e., A-B sublattice) and spin 
degrees of freedom, l ajS denotes the identity matrix in 
the a and s space, a is the Rashba SO coupling strength, 
and M is the exchange magnetization. The Hamiltonian 
for valley K' is obtained by the replacement a —> — cr*. 
We neglect intrinsic SO coupling since we are interested 
in the case when the Rashba SO coupling parameter a 
is much stronger than the intrinsic coupling parameter 
^intrinsic- We note, however, that the presence of a small 
but finite intrinsic SO coupling is not expected to qual- 
itatively modify our results as long as Ai ntr i n sic <^ M. 

Upon diagonalization of the Hamiltonian, we obtain 
the energy dispersion for both valleys 



Sfl 



H\ M 2 + e 2 k + 2a 2 + 2s J a 4 + e\ (M 2 + a 2 ), 



(2) 

where = vk, \i = ± stands for the conduction (+) and 
valence (— ) bands. Because of spin- mixing due to Rashba 
SO coupling, spin is no longer a good quantum number, 
and the resulting angular momentum eigenstates are de- 
noted by the spin chirality s = d=. The band structure 
therefore consists of two spin chiral conduction bands 



FIG. 1: (Color online) (a)-(c) Bulk states band structure, 
(a). M = 0, a = 0, (b). M/e c = 0.1, a = 0, (c). 
M/s c = 0.1, a/s c = 0.05. k c = 2it/a (a is the graphene 
lattice constant) and s c are momentum and energy cut-off of 
the Dirac model, beyond which the energy dispersion deviates 
from linearity due to trigonal warping, (d). Berry curvature 

Q — 2(Q-\ + Q__) (the factor of two arises from the two 

valleys) for a/e c = 0.05, M/s c = 0.1. ^ peaks at the k value 
where the upper valance band has its maximum and degen- 
eracy between opposite spin states occurs when SO coupling 
is absent [see panel (b)]. 



and two spin chiral valence bands. The corresponding 
eigenstates are 



Usti (k) = 

N sfl [(ie- i2 ^P sfl , ie~ 



(3) 



where <pk — tan x (k y jk x \ N s/Jj is the normalization con- 
stant 



N s , (k) = {Pi + Q% + R% + K) 2 } 



-1/2 



(4) 



and P,Q,R are functions denned as follows: 



Ps» (fc) = 



-Me\ 



a 



^ + 4 (M 2 +a 2)j (M + e fesM ) , 



. (k) = e k 
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(M + s ksfi y /2, 



R SI1 (k) = ae k (M + e ksfl ) . 



(5) 



Fig. [T] illustrates the evolution of the electronic struc- 
ture as the exchange field M and Rashba SO coupling a 
are introduced to the system. As shown in Fig. [TJb), the 
exchange field splits the original spin-degenerate Dirac 
cone into two oppositely spin-polarized copies, and this 
in turn produces spin degeneracy circles in momentum 
space at energy e = 0. Introducing the SO coupling a 
causes a gap to open up between the conduction and va- 
lence bands around this circle along which SO coupling 
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mixes up and down spins and produces an avoided band 
crossing. The momentum magnitude k = &a a t which 
the avoided crossing occurs and the gap A are given by 



MVM 2 + 2a 2 

&a = — , A 

^VM 2 + a 2 



2aM 



VM 2 + a 2 ' 



(6) 



In the insulating regime when the Fermi level lies 
within the bulk gap, the Hall conductivity a xy = Ce 2 /h 
where C is the Chern number C which can be evalu- 
ated using the Thouless-Kohmoto-Nightingale-den Nijs 
(TKNN's) formula [21]: 



(7) 



where n labels the bands below the Fermi level, and Q Tl 
is the Berry curvature 



f2 - = ^ |x| 5fc ) ' 



(8) 



with u n denoting the Bloch state for band n. Before cal- 
culating the Chern number, we briefly comment on and 
make connections with two other formulas in the litera- 
ture that are also used to calculate the Hall conductivity 
in the insulating regime. 

For two-band Hamiltonians that can be written in the 
form H = a • d, the TKNN formula can be written in the 
form 
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dd dd 

dk x dk 
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where d = d/|d| is the unit vector which specifies the di- 
rection of d. The right hand side of Eq. (|9| can be identi- 
fied as a Pontryagin index which is equal to the number 
of times the unit sphere of spinor directions is covered 
upon integrating over the Brillouin zone. For the present 
case, however, the graphene Hamiltonian contains both 
spin and pseudospin degrees of freedom, and Eq. ([9| is 
not applicable. In this case, C is given by the following 
more general form of the Pontryagin index [22] : 
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where fc M = (oj,k x , k y 



Eq. ([10]) can 
to the TKNN formula 



(10) 

), e^ v \ is the anti-symmetric tensor, 
and G = (iuj — H)~ x is the Green function. In the non- 
interacting limit we consider in this work 
be shown to be equivalent [ 
Eq. 0. 

We now evaluate the z component of the Berry curva- 
ture Q s/Jj from Eq. ([8| for the bands which are labeled by 
s and ii. For each valley, we find that the Berry curvature 
is analytically expressible in terms of an exact differential 
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FIG. 2: (Color online) Band structure at a/s c = 5 x 10 3 and 
M/s c 0.02 for (a) U/s c 0.01, (b) U/e c 0.02, and (c) 
U/s c = 0.03. (d) Band gap as a function of U at a/s c = 0.005 
and M/s c = 0.02, 0.04, 0.06 as shown in legend. 



and the Chern number per valley for each band is given 

by 
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(12) 



The upper limits of the integrand in Eq. ( 12 ) can be set 
to oo because the Berry curvature is large only close to 
valley centers. 



Computing Eq. (12), we find that in this continuum 
model the Chern number C s _ for the individual valence 
band with s = ± is not quantized but instead depends 
numerically on the specific values of a and M. We find, 
however, that the two contributions always sum to 1 and 
therefore each valley carries a unit Chern number. Tak- 
ing into account both valleys, it follows that the quan- 
tized Hall conductivity is 



2— sgn(M). 



III. GATED BILAYER GRAPHENE 
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We extend our discussion to the case of bilayer 
graphene. In the vicinity of valley K, we can write the 
bilayer graphene Hamiltonian in the presence of Rashba 
SO coupling a, exchange field M and potential imbalance 
U as (r denotes Pauli matrices for the layer degrees of 
freedom) : 



H 



vct ■ kl s + Ml a s z + ( aT + aB ) ( a x s) z 



(a x s) z + m a l s 



-t±l s (d x T x + (TyTy) , 
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(14) 



4 



M 







C = 


C = 



u 



FIG. 3: (Color online) Phase diagram of the Chern number C 
as a function of U and M. The valley Chern number occupies 
complementary regions of the phase space with C v = for 
U < M and 4sgn(£7) for U > M. 



where t± = 0.4 eV is the AB interlayer hopping energy 
[24] . For the other valley K', the Hamiltonian is given by 
the above with a — >• — cr*. For generality, we have writ- 



ten Eq. (14) allowing for different Rashba SO coupling 
strengths and for the top and bottom layers. We 
shall now set ax = «b = <^ for simplicity as the specific 
values of a on the two layers do not alter the topology of 
the bands in our discussions below. 



The Hamiltonian Eq. (14) is not diagonalizable ana- 
lytically [25]. We therefore obtain the eigenenergies and 
eigenvectors numerically and use these to compute the 
Berry curvature. Fig. [2] shows the band structure evo- 
lution for increasing values of U: U < M (panel a) , 
U = M (panel b), and U > M (panel c). For U < M, an 
inverted-gap profile appears that is similar to the single- 
layer graphene case (Fig. [TJ^). At U = M, the gap closes 
exactly at k = 0, and reopens when U > M. Fig. [2] 
shows the behavior of the gap A as a function of the po- 
tential difference U for various values of M. We find that 
A initially increases with U and then decreases toward 
zero when U approaches the value of M . after which A 
increases again with U . 

The Berry curvature Eq. ([8| can be expressed in a 
form that is more convenient for numerical computation. 
For the n th band, the Berry curvature per valley can be 
expressed as 



QJ 1 =-2 V 



lm{{n\v x \n'){n'\v y \n)} 



(15) 



where v XiV = dH/dk XyV . 



Hamiltonian Eq. 
we find that 



Numerically diagonalizing the 
14|) and computing the Chern number 



4e 2 //isgn(M), U < M 
0, U > M 



For U < M, the Chern number is twice that of the single- 
layer graphene case, corresponding to four edge modes. 
The bilayer graphene system behaves as an quantum 
anomalous Hall insulator when U < M, and exhibits 
vanishing Hall effect when U > M. The gated bilayer 
graphene system therefore has a Hall current which is 
tunable by an external gate voltage. 

The potential difference U breaks the bilayer's top- 
bottom spatial inversion symmetry. This produces a val- 
ley Hall effect in which valley-resolved electrons scatter 
to opposite sides of the sample. This can be character- 
ized by the valley Hall conductivity which is defined as 
the difference between the valley-resolved Hall conduc- 
tivities aZ„ = cr 



a!S . We find that 



4e 2 /h sgn(C/), U > M 
0, U < M 



(17) 



(16) 



Therefore, despite having a vanishing Chern number 
when U > M, the system exhibits a finite valley Hall 
conductivity 4e 2 /h. The quantum anomalous Hall and 
quantum valley Hall effects thus occupy complementary 
regions in the U — M phase space, as summarized in the 
phase diagram Fig. [3] The gated bilayer graphene sys- 
tem therefore behaves, depending on whether U or M is 
larger, either as a quantum anomalous Hall insulator or 
a quantum valley Hall insulator. 



IV. EDGE STATE PROPERTIES 

We have presented an analysis of bulk topological prop- 
erties in single-layer and bilayer graphene cases using the 
low-energy Dirac Hamiltonian. In this section, we study 
the corresponding edge state properties on a finite single- 
layer and bilayer graphene sheet, and switch to a tight- 
binding representation from which we obtain the edge 
bands numerically. The finite-size single-layer and bi- 
layer graphene sheets in our calculations are terminated 
with zigzag edges along one direction, and are infinite 
in the other direction. The Hamiltonian for single-layer 
graphene case can be expressed as 

H S LG = c\ a C ja ^it R ( S af3Xdij)-Zcl a C j(3 

(ij)a (ij)oi(3 

+ M^24 a (s z ) aa c icn (18) 

ia, 

where the first term describes hopping between nearest- 
neighbors z, j on the honeycomb lattice, the second term 
is the Rashba SO term with coupling strength £r (d^ is 
a unit vector pointing from site j to site z), and the third 
term is the exchange field M. a, j3 denote spin indices. 
The bilayer graphene case is described by the Hamil- 
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FIG. 4: (Color online) Edge state band structure of a single- 
layer graphene ribbon with M = 0.1885 and £r = 0.0471 in 
units of the near-neighbor hopping amplitude t in the tight- 
binding Hamiltonian Eq. ( |18| ). These values correspond to 
M/sc = 0.02 and a/s c = 7.5 x 1CT 3 in the low-energy Hamil- 
tonian Eq. |l]). 



tonian 



HBLG = #SLG + H SLG + t± ^2 C ia C ja 

ieT,jeB,a 

+u c l Ci « - u E c U c i°> ( 19 ) 

ieT,a jeB,a 

where #Jlg are ^ ne (T) anc ^ bottom (B) layer Hamil- 
tonians of Eq. (18), vertical hopping t± between the lay- 



ers is represented by the third term and occurs only be- 
tween Bernal stacked neighbors, and the last two terms 
describe the potential difference U applied across the bi- 
layer. The parameters in the tight-binding Hamiltoni- 
ans Eqs. (18)-(19) above are related to those in the low- 
energy Hamiltonians Eq. ([!]) and Eq. ( [l4| ) as v = 3ta/2 
(a = 1.42 A is the graphene lattice constant), a = 2£r/3, 
and M, U, t± are the same in both equations. 



Single-layer graphene case 

shows the ribbon band structures calculated from 
Inside the bulk gap, we find counter propogat- 
ing gapless edge channels at each valley that are localized 
on opposite edges of the graphene sheet. In the FigJFJ we 
plot the probability density profile of edge state wave 
functions \ip\ 2 as a function of the atom positions along 
the width of the graphene sheet for the four edge states 
labeled by A, B, C, D in Figj4] It can be seen that 
states A and C are localized along the left edge whereas 
B and D are localized along the right edge. The edge 
states labeled by A and C have the same velocity and 
propagate along the same direction along one edge [in- 
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FIG. 5: (Color online) (a)-(d): Probability density across the 
single-layer graphene sheet for the edge state wave functions 
|^| 2 of the edge states labeled A, B, C, D in Fig. [4] The inset 
is a schematic which indicates he propagation direction of the 
corresponding edge modes. 



set of FigjEJa)], whereas B and D have opposite veloc- 
ity and propagate along the other edge. The two chiral 
edge modes each carry a unit conductance e 2 /h yielding 
a quantized Hall conductivity a xy = 2e 2 /h. In the ap- 
pendix we also present an envelope function analysis of 
the edge state properties. The edge state band structure 
obtained with this approach is found to be in excellent 
agreement with the tight-binding results. 



Gated bilayer graphene case 

To study the edge state properties corresponding to the 
quantum anomalous Hall phase and the quantum valley 
Hall phase, we show the edge state band structure at a 
fixed Rashba SO coupling for the two cases M > U and 
U > M in Fig. [6j In contrast to the band structure in the 
single layer case Fig. [4j we find that the bilayer graphene 
band structure becomes asymmetric at K and K'. Within 
the bulk gap, there exists two edge bands associated with 
each valley. In Fig. [7] we show the probability density of 
the edge state wave function |^| 2 for the edge states la- 
beled from A to H in Fig.|6]for both cases. The left panel 
shows the case M > U, and we find that the edge states 
labeled by A, C, E, G are localized on one edge whereas 
B, D, F, H are localized on the other edge. This corre- 
sponds to the quantum anomalous Hall phase [Fig.[8ja)], 
where there exists four parallel chiral edge modes yielding 
a quantized Hall conductance a xy = 4e 2 /h. 

For the case U > M, the right panel of the wave func- 
tion plot Fig. [7] reveals that the four edge modes A, C, E, 
G which propagate along the same direction now become 
split between the two edges: A and C travel along one 
edge whereas E and G travel along the opposite edge. 
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FIG. 6: (Color online) Edge state band structure of bilayer graphene at fixed Rashba SO coupling strength (in units of the 
hopping amplitude t) t R = 0.0471 for (a) M > U: M = 0.1885 and U = 0.0942; (b) U > M: M = 0.0942 and U = 0.2826. 
These values correspond in the continuum Hamiltonian Eq. (14) a/s c — 7.5 x 10 -3 ; M/s c — 0.02 and U/s c — 0.01 for (a); 
M/s c = 0.01 and U/e c = 0.03 for (b). 
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FIG. 7: (Color online) (a)-(h): Probability density of the edge state wave function |^| 2 for the edge states A, B, C, D, E, F, 
G, H labeled in Fig. [6] for the case M > U (left panel) and for U > M (right panel) as a function of atom position. 



Similarly, for the edge modes that travel in the opposite 
direction, B and D propagate along one edge whereas 
E and G propagate along the other edge. This is il- 
lustrated in the schematic of Fig{8^b). Since the total 
current along one edge now adds up to zero, the Hall 
conductivity vanishes. Nevertheless, two sets of counter- 
propagating edge modes that belong to different valleys 
K and K' travel along one edge. This situation bears a 
remarkable resemblance to the quantum spin Hall effect 
where one edge consists of two counter-propagating spin 
polarized modes. Due to the broken top-bottom layer 
spatial inversion symmetry, bilayer graphene exhibits a 
quantized valley Hall effect, with a v xy = 4e 2 /h. In the 
case of single-layer graphene, such a quantized valley Hall 
effect arises when the A-B sublattice symmetry is broken; 
however there is no obvious strategy for imposing such an 
external potential experimentally. Through top and back 
gating, bilayer graphene allows for a more experimentally 



accessible way to produce the quantum valley Hall effect. 
Our present results show that the quantum valley Hall 
effect can coexist with time reversal symmetry breaking, 
provided that the breaking of spatial inversion symmetry 
wins over that of time reversal symmetry. 



IV. CONCLUSION 

In summary, we have studied the quantum anomalous 
Hall effect in single-layer and bilayer graphene systems 
with strong Rashba spin-orbit interactions due to exter- 
nally controlled inversion symmetry breaking, and strong 
exchange fields due to proximity coupling to a ferromag- 
net. For neutral single-layer graphene, we find that the 
Hall conductivity is quantized as a xy = 2e 2 /h. For bi- 
layer graphene, in which an external gate voltage can 
introduce an inversion symmetry breaking gap, we find a 
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FIG. 8: (Color online) Schematic showing the direction of 
edge mode propagation (indicated by arrows) in the (a) quan- 
tum anomalous Hall phase; (b) quantum valley Hall phase. 



quantized Hall conductivity at neutrality equal to 4e 2 jh 
when the potential difference U is smaller than the ex- 
change coupling M. This anomalous Hall effect is simi- 
lar to the quantized anomalous Hall effect [26, 27 which 
can occur spontaneously in high-quality bilayers at low- 
temperatures, but is potentially more robust because it 
relies on external exchange and spin-orbit fields rather 
than spontaneously broken symmetries. When U > M, 
the system exhibits a quantized valley Hall effect with 
valley Hall conductivity 4e 2 /h. 

Two obstacles stand in the way of realizing the quan- 
tum anomalous Hall effects discussed in this paper ex- 
perimentally. It will be necessary first of all to introduce 
a sizeable Rashba spin-orbit coupling. One possibility is 
surface deposition of heavy-nucleus magnetic atoms that 
induce large spin-orbit coupling. The exchange field that 
is also required could be introduced through a proxim- 
ity effect. From our ab initio studies [10 , an exchange 
field splitting of 56 meV and Rashba spin-orbit coupling 
2.8 meV can be obtained by depositing Fe atoms on the 
graphene surface. Another possible solution is to deposit 
graphene on a ferromagnetic insulating substrate. The 
presence of a substrate breaks spatial inversion symme- 
try, and therefore also produces Rashba spin-orbit cou- 
pling. Since the exchange field is induced through a prox- 
imity effect, layered antiferromagnetic insulators, which 
are more abundant in nature can also be used and offer 
the advantage of an enlarged pool of candidate substrate 
materials. 
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APPENDIX: ENVELOPE FUNCTION ANALYSIS 
OF EDGE MODES 

In this appendix, we present results for the edge state 
band structures using an envelope function approach 
based on the continuum Dirac model. We shall present 
only the single-layer graphene case below, as the bilayer 
case does not offer as much analytic tractibility since the 
Hamiltonian Eq. (14) is not analytically diagonalizable. 
We first calculate the envelope-function band structure 
and then compare the results with the tight-binding band 
structure. 

Consider below a graphene sheet of infinite extent in 
the x direction but finite in the y direction spanning from 
y = to y = L. From the Hamiltonian Eq. ([!]), one 
can write down the eigenvalue problem satisfied by the 



wavefunction ^ 



-£ + M V (k 

v (k x + d y ) 





>At, ^Bt> ^A|, ^B; 



e + M 
i2a 






— i2a 
-e-M 
v (k x + d y ) 



v(k x 



-dy) 
M 



* = 0. 



For zigzag-edged graphene, the following boundary con- 
ditions apply 



^At (y = L ) 4>a± (y = l) = o, 
M (v = °) = M (y = o) = o. 



(21) 
(22) 



The solution of the problem Eq. (20) admits the ansatz 
\I> = e Xy *fr. Substituting into Eq. (20), we obtain the 



energy dispersion in terms of A from the resulting deter- 
minantal equation 



e = ii {M 2 + v 2 {k 2 - A 2 ) + 2a 2 

+2 5v /« 4 + v 2 Wx - >?) ( M2 + « 2 )} 1/2 , (23) 



which in turn yields four characteristic lengths ±Ai 5 2 in 
terms of the energy e: 



Ai,2 = -xhm 

V 



e 2 + M 2 ± 2 v /e 2 M 2 + a 2 (e 2 - M 2 ) 

(24) 

Note that Ai^ in general can be complex, corresponding 
to a mixture of the edge and bulk states. The eigenvec- 
tors can be obtained as 



*(A) 



- (e + M) \v 2 (A 2 - kg) + (e - M) 2 
i2a (e 2 - M 2 ) 
i2av (s + M) (k x - A) 

-v (k x + A) \v 2 (A 2 - fc 2 ) + (e - M) 2 



(25) 

where we have left out an inessential normalization con- 
stant. The total wavefunction can therefore be repre- 
sented as a linear superposition of the constituent basis 



(20) 
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wavefunctions 
vj> = 



Ci* (Ai) e XlV + L>!* (-AO e~ XlV 
+C 2 * (A 2 ) e X2V + £> 2 * (A 2 ) e~ X2V . 



(26) 



Using the boundary conditions Eqs. (21)-(22), we obtain 
the following determinant al equation 



/(Ai)e AlL f(Xi)e~^ L f(\ 2 )e x * L /(A 2 )e-^ L 
1111 
(k x — Ai) e AlL (^ + A!)e- A ^ fe-A 2 )e A ^ (fe + A 2 ) e~ x - L 
/(Ai)(fe x + Ai) /(AO^-Ax) /(A 2 )fe + A 2 ) /(A 2 ) (fc* - A 2 ) 



0. 



(27) 



r 



where /(A) = (A 2 - A: 2 ) + (e - M) 2 . With Ai, 2 given 
by Eq. (24), the band structure e as a function of k x can 
be obtained from solving Eq. (27). We illustrate in Fig. [9] 
the resulting band structure in the vicinity of a Brillouin 
zone corner, from which it can be seen that both the bulk 
and edge bands obtained from the envelope function ap- 
proach show excellent agreement with the tight-binding 
result. *On leave from University of Texas at Austin. 
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FIG. 9: (Color online) Band structure from the envelope func- 
tion approach and tight-binding model. The Rashba spin- 
orbit strength and exchange field in the tight-binding model 
are (in units of t) £r = 0.02 and M = 0.1 in Eq. |l8]), cor- 
responding to a/e c = 2.12 x 10~ 3 and M/e c = 0.01 in the 
low-energy Dirac Hamiltonian. The width of the graphene 
sheet is L 
stant. 



119 in units of the nearest- neighbor lattice con- 
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